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Astrophysical fluids are generically turbulent and this must be taken into account for most transport processes. 
We discuss how the preexisting turbulence modifies magnetic reconnection and how magnetic reconnection 
affects the MHD turbulent cascade. We show the intrinsic interdependence and interrelation of magnetic 
turbulence and magnetic reconnection, in particular, that strong magnetic turbulence in 3D requires recon- 
nection and 3D magnetic turbulence entails fast reconnection. We follow the approach in Eyink, Lazarian 
& Vishniac 2011 to show that the expressions of fast magnetic reconnection in Lazarian & Vishniac 1999 
can be recovered if Richardson diffusion of turbulent flows is used instead of ordinary Ohmic diffusion. This 
does not revive, however, the concept of magnetic turbulent diffusion which assumes that magnetic fields can 
be mixed up in a passive way down to a very small dissipation scales. On the contrary, we are dealing the 
reconnection of dynamically important magnetic field bundles which strongly resist bending and have well 
defined mean direction weakly perturbed by turbulence. We argue that in the presence of turbulence the 
very concept of flux-freezing requires modification. The diffusion that arises from magnetic turbulence can be 
called reconnection diffusion as it based on reconnection of magnetic field lines. The reconnection diffusion 
has important implications for the continuous transport processes in magnetized plasmas and for star forma- 
tion. In addition, fast magnetic reconnection in turbulent media induces the First order Fermi acceleration of 
energetic particles, can explain solar flares and gamma ray bursts. However, the most dramatic consequence 
of these developments is the fact that the standard flux freezing concept must be radically modified in the 
presence of turbulence. 



I. PURPOSE AND OUTLINE 

The purpose of this short paper is to discuss processes 
that govern the change of the magnetic field topology in 
astrophysical fluids. We claim that it is incorrect to ig- 
nore ubiquitous astrophysical turbulence while studying 
magnetized reconnection. We also show the intrinsic and 
very deep relation between magnetic reconnection and 
turbulence. 

In what follows, we discuss the problem of astrophysi- 
cal reconnection and point out the difference that exists 
in terms of reconnection between astrophysical systems 
and their present day numerical models (§2), reveal the 
relation between magnetic turbulence and magnetic re- 
connection in §3, proceed to the discussion of the mag- 
netic turbulence and magnetic fleld wandering in §4. We 
relate the model of fast turbulent reconnection to the 
Richarson diffusion in §5 and discuss flux freezing in as- 
trophysical systems in §6. The connection with previous 
studies stressing the role of turbulence in reconnection is 
provided in §7. We discuss astrophysical implications of 
our work in §8 and the summary is presented in §9. 
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II. ASTROPHYSICAL RECONNECTION VERSUS 
NUMERICAL RECONNECTION 

It is generally believed that a magnetic fleld embed- 
ded in a highly conductive fluid preserves its topology 
for all time due to the magnetic flelds being frozen-in-'- . 
Although ionized astrophysical objects, like stars and 
galactic disks, are almost perfectly conducting, they show 
indications of changes in topology, "magnetic reconnec- 
tion" , on dynamical time scales^. Reconnection can be 
observed directly in the solar corona^, but can also be in- 
ferred from the existence of large-scale dynamo activity 
inside stellar interiors^. Solar flares^ and 7-ray bursts 
(see^iS) are usually associated with magnetic reconnec- 
tion. A lot of previous work has concentrated on show- 
ing how reconnection can be rapid in plasmas with very 
small collisional rate»iSiii, which substantially constrains 
astrophysical applications of the corresponding reconnec- 
tion models. 

A theory of magnetic reconnection is necessary to un- 
derstand whether reconnection is represented correctly 
in numerical simulations. One should keep in mind that 
reconnection is fast in computer simulations due to high 
numerical diffusivity. Therefore, if there are situations 
where magnetic flelds reconnect slowly, numerical sim- 
ulations do not adequately reproduce the astrophysical 
reality. This means that if collisionless reconnection is 
indeed the only way to make reconnection fast, then the 
numerical simulations of many astrophysical processes. 
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including those in interstellar media, which is collisional 
at the relevant scales, are in error. At the same time, it is 
not possible to conclude that reconnection must always 
be fast on the empirical grounds, as solar flares require 
periods of flux accumulation time, which correspond to 
slow reconnection. 

To understand the difference between reconnection in 
astrophysical situations and in numerical simulations, 
one should recall that the dimensionless combination that 
controls the resistive reconnection rate is the Lundquist 
number^^, defined as 5 = L^Va/^^ where is the length 
of the reconnection layer, Va is the Alfven velocity, and 
A = r]c^/4:Tr is Ohmic diffusivity. Because of the huge as- 
trophysical length-scales L^, involved, the astrophysical 
Lundquist numbers are also huge, e.g. for the ISM they 
are about 10^^, while present-day MHD simulations cor- 
respond to S" < 10"*. As the numerical efforts scale as L^, 
where is the size of the box, it is feasible neither at 
present nor in the foreseeable future to have simulations 
with realistically Lundquist numbers. 



III. TURBULENCE AND MAGNETIC RECONNECTION 

While astrophysical fiuids show a wide variety of prop- 
erties in terms of their coUisionality, degree of ioniza- 
tion, temperature etc., they share a common property, 
namely, most of the fluids are turbulent. The turbulent 
state of the fluids arises from large Reynolds numbers 
Re = LV/v, where L is the scale of the flow, V is it ve- 
locity and v is the viscosity, associated with astrophysical 
media. Note, that the large magnitude of Re is mostly 
the consequence of the large astrophysical scales L in- 
volved as well as the fact that (the field-perpendicular) 
viscosity is constrained by the presence of magnetic field. 

Observations of the interstellar medium reveal a Kol- 
mogorov spectrum of electron density fluctuations (see 
Ref. [m and Ti*) as well as steeper spectral slopes of su- 
personic velocity fluctuations (see Ref. [l^ for a review). 
Measurement of the solar wind fluctuations also reveal 
turbulence power spectrum^^). Ubiquitous non-thermal 
broadening of spectral lines as well as measures obtained 
by other techniques (see Ref [TtI) confirm that turbulence 
is present everywhere we test for its existence. As turbu- 
lence is known to change many processes, in particular 
the process of diffusion, the natural question is how it 
affects magnetic reconnection. 

To deal with strong, dynamically important magnetic 
fields Lazarian & VishniaolS [henceforth LV99] proposed 
a model of fast reconnection in the presence of sub- 
Alfvenic turbulence. It is important to stress that un- 
like laboratory controlled settings, in astrophysical situ- 
ations turbulence is preexisting, arising usually from the 
processes different from reconnection itself. We claim 
further in the paper that any modeling of astrophysical 
reconnection should account for the fact that magnetic 
reconnection takes place in the turbulent environment. 

LV99 identified stochastic wandering of the magnetic 



field-lines as the most critical property of MHD turbu- 
lence which permits fast reconnection. As we discuss 
more fully below, this line-wandering widens the out- 
fiow region and alleviates the controlling constraint of 
mass conservation. The LV99 model has been success- 
fully tested recently in Reffl9l (see also higher resolu- 
tion results in Ref. [20l ). This model is radically differ- 
ent from its predecessors which also appealed to the ef- 
fects of turbulence (see more cornparisons in section 7). 
For instance, unlike Ref. [2l| and [l^ the model does not 
appeal to changes of the microscopic properties of the 
plasma. The nearest progenitor to LV99 was the work of 
Matthaeus & Lamkin22;24^ -vvJio studied the problem nu- 
merically in 2D MHD and who suggested that magnetic 
reconnection may be fast due to a number of turbulence 
effects, e. g. m ultiple X points and turbulent EMF. How- 
ever, Ref. [23 and[2J did not realize the key role of played 
by magnetic field-line wandering, and did not obtain a 
quantitative prediction for the reconnection rate, as did 
LV99. 



IV. MODEL OF MHD TURBULENCE AND MAGNETIC 
FIELD WANDERING 

Wandering of magnetic field lines in LV99 model 
mostly depends on the Alfvenic component of magnetic 
perturbations. The exact scaling of the component, at 
least within the currently discussed turbulence models 
providing spectral indexes of —5/3 or —3/2 is not impor- 
tant, and the theoretical a possibility of fast reconnection 
was demonstrated in LV99 for turbulence with a wide 
range of spectral indexes and anisotropics. At the mo- 
ment, we feel that the theory by Goldrcich & Sridhar— 
(henceforth GS95) is the most trusted one^. Therefore 
in what follows we present our estimates using this the- 
ory. 

The GS95 theory was formulated originally for a situ- 
ation when the energy injection happened at the Alfven 
velocity va- More general relations applicable to sub- 
Alfvenic turbulence were obtained in LV99, and we use 
those below. In particular, eddies in MHD turbulence 
become anisotropic and the key relation between the the 
parallel and perpendicular scales of eddies in the Alfenic 
turbulence is given by the so-called critical balance which 
can be written in terms of parallel £\\ and perpendicular 
to magnetic field £± scales of the eddies (LV99) : 

h^L.(^^y''M-//^ (1) 

where Li is the isotropic injection scale of the turbulence 
and Ma = ul /va is the Alfven Mach number of motions 
at the injection scale. Note, that the critical balance con- 
dition is only satisfied in the system of reference oriented 
with respect to the local magnetic field, which is different 
from the usual global system of reference related to the 
mean magnetic field. Thus we avoid using wave-vectors 
in characterizing the parallel and perpendicular scales^!. 
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In terms of perpendicular to magnetic field motions 
the scaling of Alfvenic turbulence is Kolmogorov-type: 
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which may be interpreted as a Kolmogorov cascade of 
mixing motions, which are not constrained by magnetic 
fields. 

The scaling relations for Alfvenic turbulence allow us 
to calculate the rate of magnetic field spreading. A bun- 
dle of field lines confined within a region of width y at 
some particular point will spread out perpendicular to 
the mean magnetic field direction as one moves in either 
direction following the local magnetic field lines. The 
rate of field line diffusion is given approximately by 



dx 



A 



(3) 



where A,, 



£ii , ii\ is the parallel scale chosen so that 
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the corresponding vertical scale, £±,is ^ (y^)^^^, and x is 
the distance along an axis parallel to the mean magnetic 
field. Therefore, using equation ([T]) one gets 
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where we have substituted (y^)^^^ for l^. This expression 
for the diffusion coefficient will only apply when y is small 
enough for us to use the strong turbulence scaling rela- 
tions, or in other words when (y^) < (ui/i;^)^. Larger 
bundles will diffuse at a maximum rate of Li{uL/vA)'^ . 
For {y) small equation ([4]) implies that a given field line 
will wander perpendicular to the mean field line direction 
by an average amount 
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in a distance x. The fact that the rms perpendicular 
displacement grows faster than x is significant. It implies 
that if we consider a reconnection zone, a given magnetic 
flux element that wanders out of the zone has only a 
small probability of wandering back into it. This finding 
plays an important role both for magnetic reconnection 
in LV99 and heat transfer (see Ref . [28l and [29l) . 



V. RICHARDSON DIFFUSION AND LV99 MODEL 

The advantage of the classical Sweet-Parker scheme of 
reconnection is that it naturally follows from the idea 
of Ohmic diffusion. Indeed, mass conservation requires 
that the infiow of matter through the scale of the contact 
region be equal to the outflow of matter through the 
diffusion layer A, i.e. 



VA 
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FIG. 1. Sweet-Parker laminar reconnection versus LV99 tur- 
bulent reconnection. Unlike the Sweet-Parker reconnection 
model , in the LV99 model the outflow is limited by magnetic 
field wandering rather than Ohmic difl'usivity. From Lazarian 
et al.^«. 




FIG. 2. Simulations of turbulent reconnection by Kowal et 
al.— . Magnetic field lines are visualized along with the fluc- 
tuations of current density. Substantial changes of the mag- 
netic fleld directions are due to reconnection, as the driving is 
subAlfvenic and the turbulent perturbations of magnetic fleld 
are of low amplitude. 



The mean-square vertical distance that a magnetic 
field-line can diffuse by resistivity in time t is 



{y\t)) ^ At. 



(7) 



The field lines are advected out of the sides of the recon- 
nection layer of length at a velocity of order va- Thus, 
the time that the lines can spend in the resistive layer is 
the Alfven crossing time Ia = L^/va- Thus, field lines 
can only be merged that are separated by a distance 



Lx/VS, 



(8) 



where S is Lundquist number. Combining Eqs. ^ and 
^ one gets the famous Sweet-Parker reconnection rate, 

Vrec = Va/VS- 
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In LV99 magnetic field wandering determines the scale 
of the outflow A (see Figure 1). Using expressions from 
the earlier section one can obtain (LV99): 



Vrec < VA miu 



1/2 



1/2- 



Ml 
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This limit on the reconnection speed is fast, both in the 
sense that it does not depend on the resistivity, and in 
the sense that it represents a large fraction of the Alfven 
speed. To prove that Eq. ^ indeed constitutes the re- 
connection rate LV99 goes through a thorough job of 
considering all other possible bottlenecks for the recon- 
nection and shows that they provide higher reconnection 
speed. It may be that because of the complexity of the 
argument, the LV99 theory was considered with caution 
by an appreciable part of the community till the time 
when it was successfully tested in Ref. [l^ (see Figure 2). 

Below we provide a new derivation of the LV99 recon- 
nection rates which makes apparent that the LV99 model 
is a natural generalization of the laminar Sweet-Parker 
model to flows with background turbulence. The new 
argument in Eyink, Lazarian & Vishniac'^ (henceforth 
ELVll) is based on the concept of Richardson diffusion. 
It is known in hydrodynamic turbulence that the combi- 
nation of small scale diffusion and large scale shear gives 
rise to Richardson diffusion, where the mean square sep- 
aration between two particles grows as once the rms 
separation exceeds the viscous damping scale. A similar 
phenomenon occurs in MHD turbulence. In both cases 
the separation at late times is independent of the mi- 
croscopic transport coefficients. Although the plasma is 
constrained to move along magnetic field lines, the com- 
bination of turbulence and ohmic dissipation produces a 
macroscopic region of points that are " downstream" from 
the same initial volume, even in the limit of vanishing re- 
sistivity. 

Richardson diffusion (see Ref. Isil ) implies the mean 
squared separation of particles {\xi{t) — X2{t)\'^) ~ et^, 
where t is time, e is the energy cascading rate and () de- 
note an ensemble averaging. For subAlfvenic turbulence 
e « u\/{vALi) (see LV99) and therefore analogously to 
Eq. - 



one can write 



L{L/Uf'Hll 
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where it is assumed that L < Li. Combining Eqs. 
and (fTOj) one gets 



Vrec,LV99 ~ VAiL/Li^^^M 



(11) 



in the limit of L < Li. Analogous considerations allow 
to recover the LV99 expression for L > Li, which differs 
from Eq. pT|) by the change of the power 1/2 to —1/2. 
These results coincide with those given by Eq. (|9]). 

It is important to stress that Richardson diffusion ul- 
timately leads to diffusion over the entire width of large 
scale eddies once the plasma has moved the length of one 



such eddy. The precise scaling exponents for the turbu- 
lent cascade does not affect this result, and all of the 
alternative scalings considered in LV99 yield the same 
behavior. 

It is also important to emphasize that the original LV99 
argument made no essential use of averaging over turbu- 
lent ensembles. The stochastic line-wandering which is 
the essence of their argument holds in every realization 
of the flow, at each instant of time. The "spontaneous 
stochasticity" of field lines is not a statistical result in 
the usual sense of turbulence theory and does not arise 
from ensemble-averaging. The only use of ensembles in 
LV99 is to get a measure of the "typical" wandering dis- 
tance A of the field lines (in an rms sense). If one looks 
at different ensemble members of the turbulent flow, or 
at different single-time snapshots of the steady reconnec- 
tion state, then A will fluctuate. Thus, the reconnection 
rate will also fluctuate a considerable amount over en- 
semble members or over time. E.g. see Figs. 12-14 in 
Ref. [l^ But it will be "fast" in each realization and at 
each instant of time, because the mass outflow constraint 
is lifted by the large wandering of fleld lines. 

The LV99 theory, therefore, does not involve "turbu- 
lent resistivity" or "turbulent magnetic diffusivity" as 
this is usually understood. This is ordinarily meant to 
be an enhanced diffusivity experienced by the ensemble- 
averaged magnetic field (B). However, it is only if one 
assumes some scale-separation between the mean and the 
fluctuations that the effect of the fluctuations can be le- 
gitimately described as an enhanced diffusivity'^ . In re- 
alistic turbulent flows, with no scale-separation, this phe- 
nomcnological description as an effective diffusivity can 
be wildly inaccurate. LV99 makes no appeal to such con- 
cepts and, indeed, never considers the ensemble-average 
field (B) at aU. 



VI. VIOLATION OF FLUX FREEZING AND 
RECONNECTION DIFFUSION 

While the derivation of the LV99 expressions in the 
previous section may look trivial, appealing to the con- 
cept of Richardson diffusion, in reality the justification of 
the treatment is rooted in fundamental progress achieved 
recently in understanding the concept of frozen-in field- 
lines for turbulent MHD plasmas^ (see also Ref. IssI ? 

It is clear that in the presence of magnetic reconnection 
occurring densely in space the magnetic field lines can- 
not preserve their identity in turbulent flows. In fact, 
field-wandering is a process that is difficult to under- 
stand within the standard concept of fiux freezing. Note 
that line-wandering implies that every space point is a 
nexus of infinitely many distinct lines. If magnetic field 
lines behaved like elastic threads which cannot change 
their topology, turbulence would create a system of unre- 
solved magnetic knots draining the energy to small scales. 
Therefore, instead of fluid-type MHD one would get vis- 
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coelastic dynamics like Jello or rubber. 

The high speed of reconnection given by equation Q 
naturally leads to a question of self-consistency. Is it 
reasonable to take the turbulent cascade suggested in 
GS95 when field lines in adjacent eddies are capable of 
reconnecting? It turns out that in this context, our es- 
timate for Vrec. global IS just fast cnough to be interest- 
ing. We note that when considering the intersection of 
nearly parallel field lines in adjacent eddies the acceler- 
ation of plasma from the reconnection layer due to the 
pressure gradient is not i^'^^v^, but rather {i^^/£\)^^v\, 
since only the energy of the component of the magnetic 
field which is not shared is available to drive the outflow. 
On the other hand, the characteristic length contraction 
of a given field line due to reconnection between adjacent 
eddies is only {£^\/£\)~^. This gives an effective ejection 
rate of va/(^\\- Since the width of the diffusion layer over 
a length is just £±, we can replace equation ^ with 

Vrec,eddy ~ ^A^-- (12) 

The associated reconnection rate is just 

'^reconnect ~ (13) 

which in GS95 is just the nonlinear cascade rate on the 
scale Note, that this result is general and does not in- 
volve assuming that GS95 model of turbulence is correct. 
It, however, assures that the turbulent reconnection rates 
are high enough to avoid the formation of unresolved 
magnetic field knots, i.e. unresolved field line intersec- 
tions, in the GS95 and similar models of 3D turbulence. 

Another consequence of magnetic reconnection is the 
diffusion of plasma between different magnetic lines in 
turbulent flow. This process allows convective heat trans- 
fer as described in Ref. ? as well as magnetic field 
removal from molecular clouds and accretion disk o^^'^^ . 
The astrophysical diffusion enabled by reconnection of 
magnetic field lines gives rise to a concept of reconnec- 
tion diffusion, which for many processes, e.g. for various 
stages of star formation may be more efficient than the 
ambipolar diffusion. Below we briefly discuss this new 
concept (see also Ref. [13, [H, and |37|) . 

The common wisdom based on the notion of nearly 
perfect fiux freezing for astrophysical magnetic fields sug- 
gests that mixing of plasma entrained on different fiux 
tubes is impossible. However, it is easy to see that this 
is not true in the presence of LV99 reconnection. First 
of all, the process can connect magnetic field lines with 
different plasma loading, inducing diffusion and mixing 
along the newly emerging magnetic field lines. Then, in 
turbulent fluid the magnetic field lines can be shredded 
and mixed by eddies of smaller and smaller scales. These 
two processes are illustrated by separate figures below, 
but, in reality, these two processes in turbulent magne- 
tized fluid take place simultaneously. 

Figure 3 illustrates the concept of reconnection diffu- 
sion for two fiux tubes which have the same total pressure 
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FIG. 3. Reconnection diffusion: process 1. Illustration of the 
mixing of matter due to reconnection as two flux tubes of 
different magnetic field strength interact. 




FIG. 4. Reconnection diffusion: process 2. Illustration of the 
mixing of matter due to reconnection as two flux tubes of 
different magnetic field strength interact. 



Ptot = Ppiasma + Pmagn- In the absence of turbuleuce the 
two flux tubes are in pressure equilibrium and the en- 
trained plasma stays on the flux tubes. In the presence 
of 3D turbulence fiux tubes can reconnect (minimizing 
the energy of the ^-component of magnetic field) which 
allows plasma flow from one flux tube to another. This 
process is an illustration of a multi-scale process taking 
place in realistic turbulent flows. 

Figure 4 illustrates the concept of reconnection dif- 
fusion when the magnetic pressure in the flux tubes is 
the same. This is, for instance, is the case of heat ad- 
vection by turbulence. The mixing is happening as new 
magnetic flux tubes are constantly formed from magnetic 
flux tubes that belong to different eddies. In the figure 
two adjacent eddies are shown and the process is lim- 
ited to the effects of eddies of a single scale. It is clear 
that plasmas which was originally entrained over differ- 
ent flux tubes gets into contact along the new emerg- 
ing flux tubes. The process similar to the depicted one 
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takes place at different scales down to the scale of the 
smallest eddies. Molecular diffusivity then takes over. In 
the case of heat transfer small scale molecular diffusivity 
will ensure that the temperature along the newly formed 
magnetic flux tubes is the same. 

Naturally, the process illustrated by Figure 4 also hap- 
pens when the pressure of plasmas along magnetic flux 
tubes is different. LV99 process of magnetic field recon- 
nection ensures that both magnetic field and plasmas are 
diffusing at the turbulent diffusion rate. We also note 
that in the presence of forcing, e.g. gravitational forces 
acting upon the conducting gas, the diffusion will be ac- 
companied by the removal of magnetic field from the cen- 
ter of the gravitational potential due to the magnetic field 
bouyancy. It is important to understand that the process 
in Figure 4 is limited only by the velocity of the eddies. 
Therefore supersonic turbulence can induce supersonic 
mixing. 

Reconnection diffusion is due to eddies that are per- 
pendicular to the local direction of magnetic field. This 
direction, in general, does not coincide with the mean 
magnetic field direction. Therefore in the lab system of 
reference related to the mean magnetic field the diffusion 
of magnetic field and plasmas will happen both parallel 
and perpendicular to the mean magnetic field direction. 



VII. RELEVANT WORK ON RECONNECTION RATES 

Having discussed the LV99 model and new arguments 
supporting it, it is worth considering recent work on al- 
ternative approaches to calculating reconnection rates. 
Over the last decade, more traditional approaches to re- 
connection have changed considerably. At the time of 
its introduction, the models competing with LV99 were 
modifications of the single X-point collisionless reconnec- 
tion scheme first introduced by Petschek^. Those mod- 
els had point-wise localized reconnection regions which 
were stabilized via plasma effects so that the outflow 
opened up on larger scales. Such configurations would 
be difficult to realize in the presence of random forcing, 
which would be expected to collapse the reconnection 
layer. Moreover, Refs|H argued that observations of so- 
lar flares were inconsistent with single X-point reconnec- 
tion. 

In response to these objections, more recent models of 
collisionless reconnection have acquired several features 
in common with the LV99 model. In particular, they 
have moved to consideration of volume filling reconnec- 
tion, (although it is not clear how this volume filling is 
achieved in the presence of a single reconnection layer 
(see Ref. flTl)). While much of the discussion still centers 
around magnetic islands produced by reconnection, in 
three dimensions these islands are expected to evolve into 
contracting 3D loops or ropes'^- , which is broadly similar 
to what is depicted in Figure 1, at least in the sense of 
introducing stochasticity to the reconnection zone. 

The departure from the concept of laminar reconnec- 



tion and the introduction of magnetic stochasticity is 
also apparent in a number of the recent papers appeal- 
ing to the tearing mode instability to drive fast recon- 
nection (see Refs. IHl and LV99 showed that the 
linear growth of tearing modes is insufficient to obtain 
fast reconnection. More recent work is based on the idea 
that the non-linear growth of magnetic islands due to 
mergers provides large scale growth rates larger than the 
tearing mode linear growth rates on these scales. A sit- 
uation where the non-linear growth is faster than the 
linear one is rather unusual and requires further inves- 
tigation (see Ref. '43,) Since tearing modes exist even 
in a collisional fiuid, this may open another channel of 
reconnection in such fluids. This reconnection, as we dis- 
cuss below, should not be "too fast" to account for the 
observational data. 

If local turbulence is driven by release of energy from 
the magnetic field, it may result in a runaway turbulent 
reconnection process which may be relevant to some nu- 
merical simulations^!^. Alternatively, if tearing modes 
begin by driving relatively slow reconnection then a sim- 
ilar runaway might result^ 

In any case, in most astrophysical situations one has to 
deal with the pre-existing turbulence, which is the conse- 
quence of high Reynolds number of astrophysical fluids. 
Such turbulence may modify or suppress instabilities, in- 
cluding the tearing mode instability. In this paper we 
have shown that it, by itself, induces fast reconnection 
on dynamical time scales. 



VIII. ASTROPHYSICAL IMPLICATIONS 

Fast magnetic reconnection induces numerious astro- 
physical implications. In the sections above we have dis- 
cussed a new concept of reconnection diffusion that is 
likely to dominate diffusion of magnetic fields and plas- 
mas in various astrophysical environments from accre- 
tion disks to intracluster medium (see e.g. Refs. IH and 
50). The role of reconnection diffusion may be different. 
For instance, in circumstellar accretion disks and cores of 
molecular clouds it removes magnetic flux, while in more 
diffuse interstellar medium it mostly provides good mix- 
ing destroying the correlation of magnetic field strength 
with plasma density. Turbulent transport of heat, as well 
as impurities and dust is happening due to the reconnec- 
tion diffusion. 

At the same time, the LV99 magnetic reconnection 
induces a number of new effects. For instance, parti- 
cles trapped over magnetic field lines get accelerated via 
First order Fermi proces a^^i^^ . This acceleration is also 
reported in a situation when loops of magnetic field are 
formed via tearing reconnection^ The physics is the 
same, namely, in both cases magnetic fields shrink and 
the particles entrained over magnetic field lines get accel- 
erated. This provides a likely explanation of the origin of 
the anomalous cosmic rays measured by Voyager o^^'^^ as 
well as the anisotropy of cosmic rays in the direction of 
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heliotail^ reported by different groups. We would like to 
stress that the acceleration in 3D and 2D is found to be 
different^ and therefore results of 2D numerical experi- 
ments dealing with particle acceleration should be taken 
with grain of salt. 

Naturally, the acceleration of energetic particles and 
reconnection diffusion do not encompass all the possi- 
ble applications of the LV99 model of reconnection. For 
instance, Lazarian et al^ proposed a way of explaining 
gamma ray bursts appealing to turbulent reconnection. 
This idea was further elaborated in a high impact pa- 
per by Zhang & Yan^ In general, bursts and flares are 
the natural consequence of the LV99 model. If magnetic 
fields are originally laminar, the low reconnection rate al- 
lows the accummulation of magnetic flux. As turbulence 
increases, for example, due to the outflow, the reconnec- 
tion rate increases, making the outflow more turbulent. 
This induces positive feedback resulting in what can be 
termed "reconnection instability" . 

Whatever the particular astrophysical consequences of 
LV99 model, the most striking is the fact that in turbu- 
lent fluids the basic idea of flux freezing is dramatically 
modified in turbulent fluids. This shakes the foundations 
of the astrophysical MHD, as astrophysical fluids are usu- 
ally turbulent. 



IX. SUMMARY 

Our main points can be summarized as follows: 
Astrophysical fluids are turbulent and turbulence must 
be accounted in the models of astrophysical magnetic re- 
connection. 

Turbulence and magnetic reconnection are two inter- 
dependent processes: turbulence makes magnetic recon- 
nection fast, but magnetic reconnection is required for 
the turbulence to evolve in a self-similar fashion. 

Plasma effects may be important for local reconnection 
effects, but this in most cases will not affect the resulting 
global reconnection rates. 
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